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Abstract 


The Shape Calculus is a bio-inspired calculus for describing 3D 
shapes moving in a space. A shape forms a 3D process when combined 
with a behaviour. Behaviours are specified with a timed CCS-like pro- 
cess algebra using a notion of channel that models naturally binding 
sites on the surface of shapes. In this paper, the full formal timed 
operational semantics of the calculus is provided, together with exam- 
ples that illustrate the use of the calculus in a well-known biological 
scenario. Moreover, a result of well-formedness about the evolution of 
a given network of well-formed 3D processes is proved. 


1 Introduction 


The language of the Shape Calculus was previously introduced in this volume 
[2] with the aim of gently and incrementally present all its features and 
their relative semantics. The motivations behind the type and nature of the 
calculus operators were discussed and a great variety of scenarios in which 
the calculus may be used effectively were described. The Shape Calculus is 
intended to be a core language providing basic operators to describe a large 
variety of biological scenarios. A network of 3D processes, as introduced 
in [2], represents a virtual environment in which biological entities with 
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their physical 3D shape are placed. The dynamics of the calculus describes 
how they can interact through their channels, bumping or binding. In case 
of binding, the composed process can exhibit new channels and interact 
in a different way. Moreover, a composed process can split weakly - non- 
deterministically at any time on one bond - producing two pieces, or strongly 
- on a set of bonds as soon as the whole strong-split is possible - producing 
several pieces like in a biochemical reaction. 

A formal definition of the timed operational semantics of the calculus 
is needed to integrate and complete the partially informal description given 
in [2]. It is also the starting point for constructing both a suitable simulator 
for the calculus and several verification tools that can be used to prove prop- 
erties of the defined model or of an abstraction of it. In this paper, such a 
semantics is fully provided. Moreover, we present a first fundamental result 
of correctness and suitability of the calculus. A concept of well-formedness 
is introduced, starting from shapes and porting it to more complex calcu- 
lus objects such as 3D processes and, ultimately, networks of 3D processes. 
Well-formedness is a standard concept used to avoid strange and unwanted 
situations in which a term can be legally written by syntax rules, but that 
semantically corresponds to a contradictory situation. In our case, an ex- 
ample of such a situation is, for instance, a composed shape whose pieces 
move in different directions. We prove that a given well-formed network 
of 3D processes always evolves into a well-formed network of 3D processes, 
that is to say, no temporal and spatial inconsistencies are introduced by the 
dynamics of the calculus. 

The reader of this paper is strongly encouraged to read [2] first, in order 
to get familiar with the fundamental concepts of the calculus and, thus, to 
follow more easily the technical definitions and results. The proofs of several 
results are omitted here. They can be found in [1], an extended version of 
this paper that contains also a larger introduction to the general concepts of 
the Shape Calculus. The rest of the paper is organized as follows. Section 2 
introduces 3D shapes, shape composition, movement, collision detection and 
collision response. Section 3 defines behaviours and 3D processes giving 
them full semantics. Section 4 puts all the pieces together and specifies 
precisely networks of 3D processes and a general result of dynamic well- 
formedness is proven. 
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2 3D Shapes 


Let P, V = R® be the sets of positions and velocities, respectively, in a global 
three dimensional coordinate system. We also assume relative coordinate 
systems that will always be w.r.t. a certain shape S. We refer to this relative 
system as the local coordinate system of shape S. Using local coordinate 
systems we can express parts of a given shape, such as faces and vertexes, 
independently from its actual global position. If p € P is a position expressed 
in global coordinates, and V C P is a set of points expressed w.r.t. a local 
coordinate system whose origin is p, global(V,p) =V+p={x+p|xeV} 
denotes the set V w.r.t. the global coordinates. 


Definition 1 (Basic Shapes) A basic shape o is a tuple (V,m, p,v) where 
V CP is a convex polyhedron (e.g. a sphere, a cone, a cylinder, etc.)? that 
represents the set of shape points; m € Rt, p€ P andv € V are, resp., the 
mass, the centre of mass? and the velocity of o. All possible basic shapes are 
ranged over by o,0',---. We also define the boundary B(c) of o to be the 
subset of points of V that are on the surface of o*. 

Three dimensional shapes of any form can be approximated - with arbitrary 
precision - by “gluing” basic shapes on a common surface. This concept is 
generalized by the following definition. 


Definition 2 (3D shapes) The set S of the 3D shapes, ranged over by 
S,S’,--+, is generated by the grammar: S := 0 | S(X)S where o is a 
basic shape and X CP. If S = a = (V,m,p,v), we define P(S) = V, 
m(S) = m, R(S) = p, v(S) = {v} to be, resp., the set of points, the 
mass, the reference point and the velocity of S. If S = S,(X)S> is a 
compound shape, then: P(S) = P(S1) UP(S2), m(S) = m(S1) + m(S2), 
R(S) = (m(S1) : R(S1) + m(S2) . R(S2))/(m(S1) + m(S2))° and v(S) = 
v(S1) Uv(S2). Finally, the boundary of S is defined to be the set B(S) = 
(B(S;) U B(S2)) \ {x € P | x is interior of P(S1 (X) S2)}, where a point 
x € V CP is said to be interior if there exists an open ball with centre x 
which is completely contained in V. 


?From a syntactical representation point of view, we assume that. V is finitely repre- 
sented by a suitable data structure, such as a formula or a set of vertices. 

?We actually need only a reference point. Thus, any other point in V can be chosen. 

4Note that we consider only closed shapes, i.e. they contain their boundary. 

° Again for simplicity, we use the centre of mass as the reference point. Any other point 
can also be chosen. 
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In the following, we only consider shapes that are well-formed according to 
the following definition. 


Definition 3 (Well-formed shapes) Each basic shape o is well-formed. 
A compound shape S; (X) S2 is well-formed iff: 


1. both S, and S2 are well-formed; 
2. the set X = P(S,) A P(S2) is non-empty and equal to B(S,) M B(S2); 
3. v(S1(X) $2) is a singleton {v}° where v = v(S1) = v(S2). 


We say that two shapes S; and S2 interpenetrate each other if there exists 
a point x that is interior of both P(S,) and P(S2). In other terms, they 
interpenetrate iff the set X = P(S1)NP(S2) £0 is not a subset of B(S1) M7 
B(S2). If S = $1 (X) So ts well-formed, X is said to be the surface of contact 
between S,; and Sz; moreover, each x € X is a point of contact. 


Condition (2) guarantees that well-formed compound shapes touch but do 
not interpenetrate; the surface of contact X (a single point, a segment or a 
surface) is always on the boundary of both $; and Sj. Most of the time X 
is a (subset of) a feature of the basic shapes composing the 3D shape, i.e., a 
face, an edge or a vertex. Condition (3) imposes that all the shapes forming 
a compound shape have the same velocity; thus, the compound shape moves 
as a unique body. Since a compound 3D shape S can be represented in a 
number of different ways by rearranging its basic shapes and surfaces of 
contact, a structural congruence is defined to ‘equate’ all these possible 
representations. 


Definition 4 (Structural Congruence of 3D Shapes) The structural 
congruence relation over S, denoted by =g, is the smallest relation that sat- 
isfies the following rules: 


= Si (X) So =9 So (X) Si; 
- (Sy (X) S32) (Y) S3 =s S} (X) (So (Y) S3), where Y C B(S2) M B(S3). 


°With abuse of notation, throughout the paper, we write v(S) also to refer to the 
element v of the singleton {v} as this is not ambiguous when S is well-formed. 
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2.1 Trajectories of Shapes 


The general idea of the Shape Calculus is to consider a three-dimensional 
space in which several shapes reside, move and interact. Variations of shape 
velocities over time are treated by approximating a continuous trajectory of 
a shape with a polygonal chain, i.e. a piecewise linear curve in which each 
segment is the result of a movement with a constant velocity and vertices of 
the chain corresponds to velocity updates [3]. 

In detail, a global parameter A € R™, called movement time step, rep- 
resents the maximum period of time after which the velocities of all shapes 
are updated. The quantification of A depends on the desired degree of ap- 
proximation and also on other parameters connected to collision detection. 
In some situations, the time of updating can be shorter than A because, 
before that time, collisions between moving shapes can occur. These colli- 
sions must be resolved and the whole system must re-adapt itself to the new 
situation. The time domain T = R¢ is then divided into an infinite sequence 
of time steps t; such that to = 0 and t; < t;_1 +A for all i > 0. The updat- 
ing of velocities is performed by exploiting a function steer: T > (S @ V) 
that describes how the velocity of all existing shapes (i.e. all shapes that are 
currently moving in the space) at each time t is changed. We assume that, 
at any given time instant t € T, steer t S is undefined iff shape S does not 
exist and, hence, its velocity has not to be changed. 


Definition 5 (Evolution of shapes over time) Let S € S and t € T; 
S+t, ie. the shape S after t time units, is defined by induction on S: 
Basic: (V,m,p,v)+t=(V+(t-v),m,p+ (t-v),v) 

Comp: (51 (X) So) +t = (Si 4+ t) (X +t- v(S)) (So +t) 


Definition 6 (Updating shape velocity) Let S ¢ S and w € V. We 
define the shape S|w], i.e. S whose velocity is updated with w, as follows: 
Basic: (V,m,p,v)||w] = (V,m, p, w) 

Comp: (Si (X) S2)[w] = (Sillw]) (X) (Sallw]) 


The following result comes directly from Def. 3. 


Proposition 1 Let Se€S,teT andweV. If S is well-formed then S+t 
and S|[w] are well-formed. 
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2.2 Collisions: Detection and Response 


The intent to represent a lot of shapes moving simultaneously in space 
produces inevitably a scenario where collisions among shapes are possible. 
There is a rich literature on collision detection systems. For our purposes, it 
is sufficient to define an interface between our calculus and a typical collision 
detection system pivoting on the notions described in the following. 


Definition 7 (First time of contact) Let I be a non-empty finite set of 
indexes and let {S;}ier be a set of well-formed shapes such that for alli,j € I, 
5S; and S; do not interpenetrate (see Def. 3). The first time of contact of 
the shapes S;, denoted Ftoc({S;}ier), is a number t € T such that: 


1VWET.0<t' <tandVi,j € I. S;4+t' and S;+t' do not interpenetrate; 


2. Fi,jg € 1.1 #j, such that B(S; +t) B(S; +t) AO, t.e., some shapes 
are touching at t; 


3. Ve > 056.0 <6 <€ and i,j € 1.1 4 j, such that S; + (t +6) and 
Sj + (t+) interpenetrate, t.e., int some shapes are touching and any 
further movement makes them to interpenetrate. 


Such a definition allows shapes that are touching without interpenetrating, 
and with velocities that do not make them to interpenetrate (e.g., the same 
velocity), to move without triggering a first time of contact. This will be 
useful in Section 4 when we split previously compound shapes. Indeed, after 
the split these shapes have the same velocity and, hence, do not affect the 
next first time of contact. 


Definition 8 (Collision information) Let {S;}ic7 be a set of well-formed 
shapes and let t = Ftoc({S;}ier) be their first time of contact. The set of 
colliding shapes after time t is denoted by colliding({S;}ier) CSS x 9(P). 
A tuple (5;,5;,X) € colliding({Si}ier) if: 


1. P(S;+t)NP(S;+t) is non-empty and is equal to B(S;+t)NB(S; +t); 


2. Ve > 050.0 < 6 <€, such that S;+(t+6) and S;+(t+0) interpenetrate. 


The problem of collisions response [4], i.e. how collisions, once detected, 
are resolved, is treated by distinguishing between elastic collisions (those 
in which there is no loss in kinetic energy) and perfectly inelastic ones (in 
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which kinetic energy is fully dissipated)’. After an elastic collision, the two 
shapes will proceed independently to each other but their velocities will be 
changed according to the laws for conservation of linear momentum and 
kinetic energy - Equations (1)-(2) in Def. 9. On the contrary, two shapes 
that collide inelastically will bind together and will move as a unique body 
whose velocity is given by the law for conservation of linear momentum only 
- Equation (3), in Def. 9. 


Definition 9 (Collision Response) Let 51,52 € S and let X C P be a 
surface of contact. If X is neither an edge nor a vertex of S1, the velocities 
wy, and we of these shapes after an elastic collision in X are given by: 


d d 
m(Si) ee m(S2) 


(1) wi = v(S1) — 
where n is the normal of contact away from X C B(S}), i.e. the unit vector 
perpendicular to the face of S, that contains X, and 


4 m(S1)m(S2) v(S1)-n—v(S2)-n 
~~ m(S1) + m(S$2) n-n 


If X is either an edge or a vertex of S1, n is the normal of contact away from 
the shape Sz and velocities w, and W2 are obtained by means of symmetric 
equations. In both cases, we write Sy ey, S2 to denote the pair of velocities 
(wi, Wwe). If S; and S2 collide inelastically in X, they will bind together as 


a unique body whose velocity (denoted with Si ea So) is given by: 


_ mS) 
oh ais. a asaya 


m(S2) 
($1) + m(S2) 


-v(S2) 


3 3D Processes 


Now we introduce the timed process algebra whose terms describe the in- 
ternal behaviour of our 3D shapes. We briefly recall the definitions of shape 
behaviours, as well as the associated temporal and functional semantics. 


Definition 10 (Shape behaviours) The set of shape behaviours, denoted 
by B, is generated by the grammar: 


"Other different kinds of collisions can be easily added to the calculus provided that 
the corresponding collision response laws are given. 
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NIL;———— PREF; STR ; 
nil ~ nil BwB~» p.B p(L).B ~» p(L).B 
Bi Bl Bo ~ Bt >t BSB Ke B 
SUM; DEL¢ ; DEF; ; 
Bick Be Bra ee e(').B ~ €(t! —t).B Ks» B! 


Table 1: Temporal behaviour of B’s terms 


poe eet eee! 
PREFg——W\—___ DF La, SUM, 
u.B* B <(0).B 4s BY nae S 
ee L = {(a,X)} 
DEF, STR, a 
ae p(L).B "+ B 
L={(a,X)}UL' LA0 B LX), Br 
STR2 (a,X) STR3 aX) 
oe emcee CPE EB eo nay 


Table 2: Functional behaviour of B-terms 


B::=nil| (a,X).B | w(a,X).B | p(L).B | €(t).B| B+ B| K 
where (a, X) €C, L CC (non-empty) whose elements are pairwise incom- 
patible (t.e. for each pair (a, X),(8,Y) € L it is (a,X)  (8,Y)), t€ T 
and K is a process name in K. 


Definition 11 (Operational semantics of shape behaviours) The 
SOS-rules that define the temporal transition relations 6 (B x B) for 
t € T, that describe how shape behaviours evolve by letting time t pass, are 
provided in Table 1. We write B ~> B’ if (B, B’) es and BS if there is 
B' €B such that (B, B’) €~98. Rules in Table 2 define the action transition 
relations “+C (B x B) for uw € Act, describing which basic actions a shape 
behaviour can perform. 


Now we are ready to define 3D processes, i.e. simple or compound 
shapes whose behaviour is expressed by a process in 


*Similar conventions will apply later on. 
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Definition 12 (3D PEocesses) se i. 3DP of 3D processes is generated 
by the following grammar: P : B] | P(a,X)P whereS €S, BEB, 
aéA and X is a non-empty a of P. The shape of each P € 3DP is 
defined by induction on P as follows: 

Basic: shape(S[B]) = S$ 

Comp: shape(P (a, X) Q) = shape(P) (X) shape(Q) 


We also define v(P) = v(shape(P)) and B(P) = B(shape(P)). Below we of- 
ten write P is, Q and P Poa Q as shorthand for shape(P) is, shape(Q) 
and shape(P) ree shape(Q), resp. Finally, P\v] is the 8D process we ob- 
tain by updating P’s velocity as follows: 

Basic: (S[B])|[v] = (S[v])[] 

Comp: (P(a,X) Q)|v] = (Ply]) (2, X) (QI[fv]) 


We finally write steer t P to denote Plsteer t shape(P)]|. We say that a 
basic process S|B] is well-formed iff the shape S' is well-formed and, for each 
X C P that occurs in B, global(X,R(S)) C B(S). A compound process 
P (a,X)Q is well-formed iff P and Q are well-formed, v(P) = v(Q) and 
the site X expressed w.r.t. a global coordinate system is a non-empty subset 
of B(P)NB(Q). This also means that P(P)AP(Q) is non-empty and equal 
to B(P) 1 B(Q). 


Later on we only consider well-formed processes. The following proposition 
is an easy consequence of well-formedness of shapes and 3D processes. 


Proposition 2 For each P € 3DP well-formed, shape(P) is well-formed. 


Let us recall a running example, taken from [2], based on the well- 
known biochemical pathway of glycolysis. We consider only one reaction, 
catalyzed by the Hexokinase enzyme (HEX). GLC, G6P, ATP and ADP are 
metabolites. The transformation of a metabolite (GLC) into another (G6P) 
depends on the meeting (collision in binding sites) of the right enzyme (HEX) 
with the right metabolites (GLC and ATP). After this binding the reaction 
takes place and the products (G6P and ADP) are released. A special case 
occurs when the enzyme has bound one metabolite and an environmental 
event makes it release the metabolite and not proceed to the completion of 
the reaction. We denote by S;,, Sg and S, the shape of HEX, GLC (and 
G6P) and ATP (and ADP), respectively. 


Al 


Example 1 (3D Processes for HEX, GLC and ATP) An Hexokinase 
molecule can be modelled as HEX = S;,[HEX] where HEX = (atp, Xpq)-HA+ 
(glc, Xhg)-HG, HA = w(atp, Xpq). HEX + e(tn).(gle, Xng). 

p({(atp, Xna), (gle, Yang) }). HEX, HG = w(gle, Xp,).HEX + e(t;,).(atp, Xna)- 
p({(atp, Xna), (gle, Yng) })-HEX, and Xha, Yng are proper surfaces of contact. 

ATP =S,[ATP] models an ATP molecule where ATP = (atp, Xap). 
(€(ta).p({ (atp, Xan) })-ADP + w(atp, Xa,).ATP) and the surface of contact 
Xan is the sheles boundary B(S,). The process modelling a molecule of glu- 
cose is similar: GLC = S,{[GLC] where GLC = (glc, Xgn).(€(tg). 
p({(glc, Xgn) }).-GOP + w(gle, Xgn).GLC. We leave unspecified the behaviours 
GOP and ADP. 

HEX has two binding capabilities along the channels (atp, Xha) and 
(glc, Yng) to bind, resp., with an ATP and a GLC molecule. By performing 
an action (atp, Xna), HEX evolves in HA. HA can perform either a weak- 
split action w(atp, Xpq) to come back to HEX, or can wait at most tp, units of 
time, perform (glc,Yng) and then evolve in p({(atp, Xna), (gle, Yng) }).HEX. 
Now, two strong-split actions are enabled after which we come back to HEX. 
Note that, after an action (glc, Yng), HEX becomes HG that behaves symmet- 
rically. An ATP molecule performs a (atp, Xqn)-action, waits t, units of 
time, and then can release the bond established on the channel (atp, Xan) — 
and thus return free as ATP — or can participate in the reaction and become 
an ADP. As we will see in Section 4, the result is the split of the complex 
in the three original shapes whose behaviours are HEX, ADP and G6P, resp. 
We omit the description of the behaviour of GLC since it is similar to that 
of ATP. 


The timed operational semantics of 3D processes is defined below. 


Definition 13 (Transitional semantics of 3D processes) Rules in Ta- 
ble 3 define the transition relations ae (3DP x 3DP) for t € T, and 
“,¢ 3DP x 3DP for pw € Act. Two 3D eed and Q are said to 
be compatible, written P ~ Q, if P ao and Q scan for some compatible 
channels (a, X) and (@,Y); otherwise, P and Q are incompatible, written 
PQ. We often write P dy and P 7+ as shorthand for P De ar 

P oA Ny resp., for any (a, X). 

The following proposition shows that 3D processes well-formedness is closed 
w.r.t. transitions ~> and 4 
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BY B’ Pw P!QwQAX'=X+(t-v(P)) 


BASIC; ; COMP; ; 
S[B] ~ (S +t)[B] P (a, X)Q ~ P' (a, X') Q! 
B&, BI y = global(X,R(S)) pee 
BASIC, (a¥) COMP, (n¥) 
S[B] ——> S[B" P (a,X)Q = P’ (a, X)Q 


B LO), BI Y = global(X,R(S)) 


BASIC, 


s[B] 22, s[B7 
pais p &), py C B(P (a, X) Q) 
COMP, wa¥) COMP, Ta¥) 
P (a, X) Q ——> P' (a, X)Q P (a, X) Q ——> P’ (a, X)Q 


Table 3: Functional and temporal behaviour of 3DP-terms 


Proposition 3 Let P € 3DP well-formed. Either P Re Q or P uy Q implies 
Q € 3DP well-formed. 


Notice that the operational rules in Table 3 do not allow synchroniza- 
tion between components of compound process that proceed independently 
to each other. Consider, as an example, P = S,[p({a,Xp}).Bp], Q = 
Sqle({@, Xq}).Ba], and P (a,X)Q where X = Xj, Xj and, for each i € 
{p,q}, Xi = global(X;,R(S;)), ie. Xf is the site X; w.r.t. the global co- 
ordinate system. As stand-alone processes, P and Q can perform two 
compatible strong-split actions, namely p(a,X;,) and p(@, Xj) and evolve, 
resp., in S,[B,] and S,[B,]. As a consequence, P (a, X)Q becomes either 
Sp[Bp] (a, X) Q or P (a, X) Sp[Ba]. 

Being these actions compatible, P and Q have to synchronize on their 
execution in order to split the bond (a, X) (i.e., a strong-split operation 
is enabled). Such an operation must be performed before time can pass 
further and must produce as a result two independent 3D processes, i.e. 
the network of 3D processes (see Section 4) that contains both S,[B,] and 
S,[B,]. Similarly, we would allow synchronizations between compatible 
weak-split action in order to perform a weak-split operation. To properly 
deal with this kind of behaviours some technical details are still needed. 


43 


We first allow synchronization on compatible split actions by introduc- 


i ae ; x x da 
ing the transition relations mas and gis : Intuitively, we want that 


P (a, X)Q ma) S»|Bp] (a, X) Sg[Bg]. Now, we can ‘physically’ remove the 
bond (a, X) (this will be done by exploiting the function split over 3D pro- 


cesses we provide in the next section) and obtain the network of processes 
we are interested in. 


Definition 14 (Semantics of strong and weak splits) The SOS-rules 


that define the transition relations Hee 3DP x 3DP where p(a, X) € p(C) 
are given in Table 4. We omit symmetric rules and those ones defining the 
transition relations MES) G 3DP x 3DP for w(a, X) € w(C), since these can 


be obtained from those in Table 4 by replacing each action p(-) with the 
corresponding action w(-). 


Hp OPE Ot ae eat. Se 
STRSYNC (aX) 
Pig, X) QoS" Pia, X).Q! 
pris) pi 
STRPAR (b,Y) 
P(a,X)Q°S P’ (a, X)Q 


Table 4: Transitional semantics for strong-split actions 


Recall that strong-split operations require simultaneous split of multiple 
bonds. In this case, all the components involved in the reaction must all 
together be ready to synchronize on a proper set of compatible strong-split 
actions. This concept is formalized by the definition below. 


Definition 15 (Bonds of 3DP-terms) The function bonds : 3DP > ¢(C) 
returns the set of bonds that are currently established in P. It 1s defined by 
induction on P € 3DP: 

Basic:  bonds(.$[B]) = 0 

Comp: bonds(P (a, X) Q) = bonds(P) U bonds(Q) U {(a, X) } 


Xx) . : 
By induction on P we can prove that P es implies (a, X) € bonds(P). 
Moreover, we say that P € 3DP is able to complete a reaction, written P \, 


iff either (1) P #, or (2) P aes) Q for some p(a, X) and Q such that Q \. 
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Finally, if P is able to complete a reaction and there is at a least a bond 
that has to be strongly split (i.e. if P # does not hold), a reaction can 
actually take place and, as a consequence, time cannot pass further. Below 
we restrict the timed operational semantics of 3D processes as it has been 
defined in Def. 13 in order to take this aspect into account. 


Definition 16 Let P € 3DP. We say that P Bt Q if P Be Q and either 
P 4 or P is not able to complete a reaction. 


Proposition 4 states that the function split is well-defined up to structural 
congruence over 3D processes we define below. 


Definition 17 (Structural congruence over 3D processes) We defi- 
ne the structural congruence over processes in 3DP, which we denote by =p, 
as the smallest relation that satisfies the following axioms: 


- S[B] =p S'[B] provided that S = S'; 

- P (a, X)Q =p Q (a, X) P; 

- P=pQ implies P (a, X) R=p Q (a, X) R; 

- Y C B(Q)NB(R) implies (P (a, X) Q) (b,Y) R=p P (a, X) (Q (b,Y) R). 


Proposition 4 Let P € 3DP well-formed. If (a, X) € bonds(P) there is a 
well-formed 3D process Q (a, X) R=p P. 


We also need the following closure result. 


Proposition 5 Let P,Q € 3DP and pw € w(C)U p(C). Then: 
1. PSQ implies shape(Q) = shape(P); 
2. P well-formed and P a Q implies Q well-formed. 


4 Networks or 3D processes 


In the following, we recall the definition of a 3D network, 1.e. a collection of 
3D processes moving in the same 3D space. 


Definition 18 (Networks of 3D processes) The set N of networks of 
3D processes (8D networks, for short) is generated by the grammar: 

N ::= Nil | P | N||N where P € 3DP. Given a finite set of indexes I, 
we often write (|| P;)ier to denote the network that consists of all P; with 
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i € I. We assume that I = 0 implies (|| P;)icr = Nil. For N = (|| Pi)ier 
we let S; = shape(P;), fori € I, and define colliding(NV) as the set of all 
tuples (P;,P;,X) such that (S;,5;,X) € colliding({Si}ier) (see Def. 8). A 
network N = (|| Pi)ier is said to be well-formed iff each P; is well-formed 
and, for each pair of distinct processes P; and P;, the shapes S; and S; do 
not interpenetrate. Moreover, steer t (|| P;)ier = (|| steer t P;)ier. 


Definition 19 (Splitting bonds) The function split : 3DP x e(C) > N 
is defined as follows: If (a,X) € bonds(P) NC and P =p Q(a,X) R then 
split(P, C) = split(Q, C) || split(R, C); if, otherwise, bonds(P) 1 C = 0, then 
split(P, C) = P. 


It is worth noting that split shapes maintain the same velocity until the next 
occurrence of a movement time step. As we mentioned above, this is not a 
problem because they will not trigger a collision and, thus, a shorter first 
time of contact. 


Proposition 6 Let P € 3DP well-formed and C CC. Then split(P,C) is a 
well-formed network of 3D processes. 


Definition 20 (Semantics of weak- and strong-split operation) Let 
P € 3DP a 8D process. If P \, we write that P *NeEN iff there 
is a non empty set C = {(a1,X1),--+ , (Gn, Xn)} C bonds(P) such that 


P= Pp Bisa RO pe pe wd N= Spl O)i. Se 


ilarly, P =>» N iff d(a,X) € bonds(P) such that P a) Q and N = 
split(Q, {(a, X)}). 


Since weak-split operations are due to a synchronization between just a pair 
of 3D processes, condition ‘P is able to complete a reaction’ is not needed, 


x 
but ‘P pec ) Q’ suffices to our aim. 


Example 2 Let us consider the 8D process P = H (atp, X) (A(glc, Y) G) 
where H = Sp[ p({(atp, Xna), (gle, Yng) })-HEX], A = Sal p({(atp, Xan) })- 
ADP + w(atp, Xan)-ATP],G = S,[p({ (gle, Yon) }).G6OP + w(glc, Ygn).GLC], 
Xia 1X1, = X (here X}, and X!, are the sites Xnq and Xap, expressed 
w.r.t. a global coordinate system; this convention will be applied later on) 
and Y;, Yj, =Y. According to the definitions given so far, P is able to 
complete a reaction since: 


P PCB" Salo({(ele, Yag)})-HEX] (atp, X) (Sal[ADP] (gle, Y) G) 


nes), [HEX] (atp, X) (Sa[ADP] (gle, Y) Sp[G6P]) = R 
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Moreover, R LS and 


split(R, C) Sa {HEX] || split(S,[ADP] (gle, Y) Sg[G6P], C) 


Sp [HEX] || (Sa[ADP] || Sg[G6P]) = N, 


where C = {(atp, X), (glc, Y)}, implies P ’..N. Moreover, let us note that, 


for eacht ET, P ~~» P but P ai since P is able to complete a reaction and 
P is does not hold. 


Below we recall the temporal and functional behaviour of 3D networks. Sym- 
metric rules have been omitted. 
Definition 21 (Temporal and Functional Behaviour of N-terms) 
The following rules defines the transition relations a NxN fort € T and 
Vv 
>CNxXN forv € {w, p}: 
N->N' M-+M’ N + NI 
EMPTY;—————___ PAR; PARg 
Nil “> Nil N || M +> N’|| M’ N|| M4 N’||M 


A timed trace from N is a finite sequence of steps of the form N = No —> 
Ni--- “> Ny, = M. We also write that N aN M if there is a timed trace 


n 
N=No— N,--- => Ny =M such that t = > {v;|% € T}. 
i=0 


Proposition 7 Lett ¢T, PE 3DP, NEN, with P and N well-formed. 
1.P%N andP4N implies NEN well-formed. 
29NS+M implies M well-formed; 


3NSM implies M well-formed. 


4.1 Collision response 


In this section we describe the semantics of collisions response. The notion of 
compatibility between channels (and, hence, processes) has been introduced 


to distinguish between elastic and inelastic collision. In particular, collisions 
xX 
among compatible processes are always inelastic. So, if P ase P’ and 


a,x ’ = ‘ Hed 

Q aU Q’, with (a, X,) and (a, X,) compatible, and P and Q collide in the 
non-empty site X = X,MXq we get a compound process (P’ (a, X) Q’) |v] 
where the velocity v is provided by Equation (3) in Def. 9. Vice versa, a 


collision between two incompatible processes P and Q is treated as an elastic 
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one. After such a collision, P and Q (actually the processes we obtain by 
updating their velocities according to Equations (1) and (2) in Def. 9) will 


proceed independently to each other. To resolve collisions, we introduce two 


P,Q,X 
different kinds of reduction relations over 3D networks, namely RDAs) and 


P,Q,.X : 
eee) i, where P,Q are 3D processes and X is a surface of contact (see 


Table 5). Intuitively, if N 22%, a (w 2), Mr), then M is the 3D 
network we obtain once an elastic (inelastic, resp.) collision between P and 
Q@ in the surface of contact X has been resolved. These reduction relations 
also use the structural congruence over 3D networks. 


Definition 22 (Structural congruence over 3D networks) 
The structural congruence over terms in N, that we denote with =, is the 
smallest relation that satisfies the following axioms: 


- Nill|| N=N,N||M=M||N and N || (M|| RB) = (N || M) || B; 
- P||N=Q||N provided that P =p Q. 


Rule ELAS in Table 5 simply changes velocities of two colliding but incom- 
patible processes guided by Equations (1) and (2) in Def. 9, while rule INEL 
joins two compatible processes P and Q to obtain a compound process whose 
velocity is given by Equation (3) in Def. 9. We force P and Q to synchronize 
on the execution of two compatible actions (a, Xp) and (@, X,) before joining 
them. In rules ELAS= and INELAS= we also consider structural congruence 
over nets of processes. In Def. 23 we collect together all the reduction-steps 
needed to solve collisions listed in a given set of collisions colliding(V); clearly 
N is a generic 3D network. 


Definition 23 (Resolving collisions) Let N € N and (P,Q,X) a tuple 


in colliding(N). N “220, ar i either P~ Q and N 84, M or 
(P,Q,X) 


P &Q and N —*=5,. M. Moreover, we write that N + M if either 
colliding(N) = 0 and N = M or colliding(N) 4 @ and there is a finite 


P1,Q1,X PriQh,X 
sequence of reduction steps N = No ee Ny::: Reed ire Ny = M 


such, that: 
1. (P;,Q;, Xi) € colliding(Nj_1) for each i € [1, k]; 


2. colliding( Nz) = 0. 
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xX 
P->.Q= (Vp, Vq) 


(P||Q)||N 22... (Plvp] || Qlival) || V 


ELAS 


(a, Xp) (a,Xq) Xpnx 
————— oF 


P' Q Q’ aef{aa} P 3"; Q=v 


(P||Q) | N eke), (Pr (a, Xp Xq) QV) || N 


P 
INEL 


(P,Q,X) , {P,Q,X) 


N=N N———5.M N=WN' N ,M 
ELAS= INEL= 
N (P,Q,X) .M N (P,Q,X) ,M 


Table 5: Reaction rules for elastic and inelastic collisions 


Let also note that, at any given time ft, colliding(V) can be obtained from 
the set of all the pairs of processes in N that are touching at that time. 
This set and hence colliding(V) is surely finite and changes only when we 
resolve some inelastic collision (this is because, after an inelastic collision one 
or more binding sites can possibly become internal points of a compound 
process, and hence are not available any more). Moreover a collision between 
pairs of processes with the same shape can not be resolved twice. This is 
either because two processes P and Q have been bound in a compound 
process as a consequence of an inelastic collision, or because P and Q collide 
elastically and their velocities have been changed according to Equations (1) 
and (2) in Def. 9 in order to obtain two processes that do not collide any 
more. Thus, we can always decide if there is a finite sequence of reduction 
steps that allows us to resolve all collisions listed in colliding(V) and hence 
obtain a network M with colliding(M) = 0. 


Proposition 8 Let N € N, P,Q € 3DP and X a non-empty subset of P. 
Then N well-formed and N eee M implies M well-formed. 


Iterating Proposition 8 it is also: 
Lemma 1 Let N €N well-formed and N > M. Then M is well-formed. 


We are now ready to define how a network of 3D processes evolves by per- 
forming an infinite number of movement time steps. 
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Definition 24 (System evolution) Let N,M € N andt,t/ € T. We say 
that (N,t) 4 (M,t+t’) iff one of the following conditions holds: 
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1. t' =Ftoc(N) < A and N 4. in’ 5 N"” and M = steer (t+t’) N”; 
2. t =A <Ftoc(N) and NS N’ and M =steer (t +t’) N’. 


A system evolution is any infinite sequence of time steps of the form: 
, ; i-1 : i ‘ 
(No, 0) = (Nj, t1) 3 (No, ti + ta) on (Ni-1, >> ty) = (Ni, St) Ean oon 
j=l j=l 


Note that, for each 7 > 1, t; = min{Ftoc(N;_1), A}. Moreover, in order to 
make sure that processes will never interpenetrate during a system evolution, 
if Ftoc(N;-1) < A, we first resolve all the collisions that happen after time 
t; = Ftoc(N;-1) (by means of transition “+) and then apply the changes 
suggested by the function steer as described in Section 2.1. 


Example 3 This example shows a possible evolution of the 8D network 
No = (HEX || ATP) || GLC where HEX, ATP and GLC are the 8D pro- 
cesses of Example 1. Below we use the following notation: 


- H(t) = (S_,+t)[HEX], A(t) = (Sa +t)[ATP] and G(t) = (S, + t)[GLC] 
for eacht € T. Note that HEX = H(0), ATP = A(0) and GLC = 
G(0); 

- C = p({(atp, Xna), (gle, Xng)})-HEX and, for any t < tn, HA(t) = 
w(atp, Xpq). HEX + e(tp — t).(gle, Xng).C; 


- AH(t) = w(atp, Xan).ATP + e(ta — t).p({atp, Xan}).ADP for any t € T 
with t < tg; 


- GH(t) = w(gle, Xgn).GLC + e(t, — t).p({glc, Xgn}).G6P for any t < ty. 


Let t; = Ftoc(No) and assume ty < A. By the operational rules, it is 
No 4 H(ty) || A(t) || G(ti) = Ni. We also assume that colliding(Ni,) = 
{(H(t1), A(t1), X)} where X = X},,90X', #0, Xi, = global(Xna, R(S;, + 
t1)) and X!, = global(Xan,R(Sa t+ ti)). Then: Ni % P(t1) ||G(ti) = NO 
where P(t,) = ((Sp + t1)[HA(O)](atp, X) (Sa + t1)[AH(0)]) [Vaal] and Via = 
AG yes A) 
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Finally: — (No, 0) 4 (Ni,t1) where N; = steer t, Nj = steert; P(t1) || 
steer ty G(t1) = P(t1)|[vi] || G(t1) [va]. Note that: 


P(ti)lvil = ((Sh + t1)[HA(O)](atp, X) (Sa + t1)[AH(0)]) [vil 
= (5j,[HA(0)](atp, X).5;[AH(0)]) 


where Sp = ((Sp+ti)|[vi]) and St = ((Sa+ti)|[vi]). Moreover, G(t1) [v2] = 
((Sg + t1) [v2] [GLC] = $5 [GLC]. 

Let tz = Ftoc(N,) and assume to = th < min{tg, A}°. Below we write 
G' (tz) and P'(t2) to denote, respectively, the 3D processes (55 +t2)[GLC] and 
((S},+t2)[HA(tn)] (atp, X +t2-v1)(S1+t2)[AH(t2)]). Again by the operational 
rules, Ny 3 P'(t) || G'(t2) = Nf. 

Let colliding( Nj) = {(P"(t2),G'(t2),Y)} where Y = X,,9 Xj, # 
0, Xhg = global(Xng,R(S;, + t2)) and Xon = global(Xg,,R(Sj + t2)) C 
B(P'(t2)). If Von = P' (to) 4; G'(t2), then 
Ni +s (P(t2) (gic, Y) G(t 2)) [vgn] = Ni where G(t2) = (Sj + t2)[GH(0 )] and 
P(tz) = ((S;, + te)[C(atp, X + ta - vi) (Sq + t2)[AH(t2))). 
Finally: —(Ni,t,) 3 =, (No,ti + to) where No = steer (t; + te) Ny! = 
(P(t2) (gle, Y) G(t2)) [v3]. Observe that: 
(P(ts) (ale, ¥) G(t2)) [vs] = (P(te) Evs]) (gle, ¥) (G(t2)kvs]) = 
((S2{C] (atp, X + ta - v1) S2[AH(ta))) (gle, Y) $2[GH(0)) 
where Sj, = (Sj, + ta) [va], $2 = (Sa + ta)[lvs] ‘and S53 = (55 + t2)|[vs]. 
At this stage the network contains just one process “and: as a consequence, 
no collisions are possible. Thus, Ftoc(N2) = co. Assume tg = ta—te < <i 
If we let Xq = (X +to-vi)+tg- v3 and Y, =Y +t, - v3, then: 
Na * ((S% + ty)[C] (ate. Xp) (S2 + ty)IAH(ta)) (ele, ¥o)(S3 + ty)IGH(ty)] & 
((Sj + tg)[HEX] || (Sz + t,)[ADP]) || (Sg + ty) [G6P] = 
((Sj, + A)[HEX] || (SZ + A)[ADP]) || (Sj + A)[G6P] = 
Thus: (No, ti + ta) & (N3, ti +tg+ A) where N3 = steer (ty +to+ 
A) Nj = (S#[HEX]|| $3(ADP)) || $2(G6P], S? = (52 + A)[vi] and vi = 
steer (t) +t 2+ A) (S2+A) for each i € {h,a,g}. 


°If were t2 < tp. the 3D processes HA(tz) and GLC would be no more compatible, and 
a collision between them would be treated as elastic. On the other hand, if were tz = ta 
the idling time for AH(ta) would be over. As a consequence, time would pass further only 
after the execution of a weak-split operation that splits the bond between the Hexokinases 
and the Atp molecules 

l0Tf were tg A ta—te the reaction could never proceed since the involved molecules would 
never be able to release - all together - the bonds. Thus the system would deadlock. 
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We can prove the following basic property stating that any system evo- 
lution does not introduce space inconsistencies like interpenetration of 3D 
processes or not well-formed processes. 


Theorem 1 (Closure w.r.t. well-formedness) Let N € N well-formed. 
If (N,t) 4 (M,t+t’) then M is well-formed. 


K 


Proof: Assume that (N,t) 4 (M,t +t’) because of N 4 nt SS NY and 
M = steer (t+ t’) N” (the other case — see Def. 24 — is similar). Then, by 
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Proposition 7 (item 3) and Lemma 1, N well-formed and N & N’ SN” 
implies N” and hence M = steer (t + t’) N” well-formed. 


5 Conclusions and Future Work 


We have formally defined the full timed operational semantics of the Shape 
Calculus. Moreover, a result of well-formedness has been proven: any well- 
formed 3D network will evolve in a well-formed 3D network. Several tech- 
nical proofs of the results are omitted in this paper, but they can be found 
in the extended version [1]. On top of this formal framework, we intend to 
define abstractions of the semantics in order to be able to prove, statically, 
properties of a given 3D network. Moreover, we intend to add information 
to the basic calculus, for instance by equipping any process with a steer 
component, rather than having an abstract function for movements. 
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